Spreading and localization of wavepackets in disordered wires in a magnetic field 
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We study the diffusive and localization properties of wavepackets in disordered wires in a magnetic 
field. In contrast to a recent supersymmetry approach our numerical results show that the decay 
rate of the steady state changes smoothly at the crossover from preserved to broken time-reversal 
symmetry. Scaling and fluctuation properties are also analyzed and a formula, which was derived 
analytically only in the pure symmetry cases is shown to describe also the steady state wavefunction 
at the crossover regime. Finally, we present a scaling for the variance of the packet which shows 
again a smooth transition due to the magnetic field. 



In recent years, considerable progress has been made 
in understanding the structure of eigenstates of quasi- 
one-dimensional (quasi-lD) disordered systems (see ||] 
and references therein). Most of the investigations have 
been concerned with the cases of completely preserved 
and totally broken time-reversal symmetry. The tran- 
sition between these limits, however, has been studied 
much less and little is known about it. The only results 
for the transition were based on a heuristic approach 
by Bouchaud ||, a semiclassical analysis by Imry and 
Lerner || and numerical studies based on transfer ma- 
trices by Pichard et al. Q. The main prediction of the 
above studies is the doubling of the localization length in 
the limit of strong magnetic fields while a smooth transi- 
tion towards this asymptotic limit was assumed. Similar 
results were obtained for quantum chaotic models |§. 
Recently, the doubling of the localization length was ob- 
served in sub- micron thin wires of doped GaAs 0] , where 
for increasing magnetic field strengths a continously de- 
creasing activation energy was reported, which saturated 
indeed at half of its field free value. Motivated by this ex- 
periment, the first attempt to study such a transition us- 
ing the supersymmetry technique was undertaken in . 
These authors predict that the decay of wavefunctions in 
disordered wires in weak magnetic fields, is character- 
ized by two localization lengths; the far tails decay with 
the length Z ta iis characteristic for completely broken time- 
reversal symmetry, while at shorter distances the decay 
length is ; orig = Ztails/2. 

The predictions of Jt],|| imply two different tempera- 
ture regimes in the hopping conductivity separated by a 
boundary which depends on the magnetic field strength. 
A direct numerical test of the above prediction by means 
of the transmission approach was carried out recently 
in ||. The results are contradicting the two-scale be- 
haviour of the wavefunction in the crossover regime. The 
following possibilities were proposed in order to explain 
this discrepancy pHI: (a) the Borland conjecture (which 
connects the asymptotic decay length of the transmit- 
tance to the asymptotic decay length of the wavefunc- 



tion) breaks down in the crossover regime between the 
preserved and broken time reversal symmetry, and (b) 
the two-scale localization phenomenon is due to anoma- 
lously localized states that are irrelevant for a typical 
wire. The latter case appears more likely since in 
the average wavefunction was calculated whereas, in £ ] 
under the assumption of the Borland conjecture the au- 
thors calculated the average of the logarithm of the wave- 
function. We are going to show that none of the above 
reasons can resolve the contradictions! 

In this Letter the problem of the crossover behaviour, 
from preserved to broken time reversal symmetry, is ad- 
dressed from a totally different perspective. In particular 
we study the evolution of initially (S-like wavepackets in 
a quasi-lD geometry under the influence of a magnetic 
field. We then directly calculate the logarithmic and the 
arithmetic average of the steady state wavefunction. We 
show that the extracted decay rates change smoothly as a 
function of the magnetic field for both observables. More- 
over, we study the whole steady state distribution and 
its fluctuations as a function of the magnetic field. From 
our numerical data we extract an analytical formula for 
the asymptotic wavefunction in the crossover regime. Fi- 
nally, we present arguments according to which the dif- 
fusive constant scales smoothly with respect to the mag- 
netic field. This imposes a smooth scaling behaviour of 
the variance of the packet. Our predictions are confirmed 
by extensive numerical calculations. 

The mathematical model we consider is the time- 
dependent Schrodinger equation on a ID lattice, 



dc n (t) 
dt 
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where c n (i) is the probability amplitude for an electron 
to be at site n and H nm — H® m + iaA nm is a complex 
Hermitian Band Random Matrix (BRM), which is de- 
composed into a real symmetric matrix H° and a real 
antisymmetric matrix A with imaginary weight ia. The 
entries of the two matrices are independent Gaussian ran- 
dom numbers with variance a 2 = 1 + S nm (where 5 nm 



1 



is the Kronecker symbol) if \n — m\ < b and zero other- 
wise. The parameter b defines the hopping range between 
neighbouring sites, or, in the quasi- ID interpretation, the 
number of transverse channels along a thin wire. To re- 
late the parameter a to the magnetic field B, we note 
that the perturbation of the levels is proportional to the 
magnetic flux <f> = kBbloo through an area limited by the 
localization length. Here A: is a dimensionless constant of 
order unity that depends on the specific geometry of the 
disordered wire p 10 1. Hence we expect a ~ $/®o, where 
$0 = h/e is the elementary flux quantum. 

We have integrated Eq. (Q) numerically using a Cayley 
scheme in order to preserve the norm jLl"| . Moreover a 
self-expanding algorithm was implemented to eliminate 
finite-size effects. Whenever the probability of finding 
the particle at the edges of the lattice exceeded 10 -15 , 
20b new sites were added to each end. Since all eigen- 
states of the Hamiltonian H nm are known to be exponen- 
tially localized with a localization length (E) ~ b 2 [[l3) , 
the evolution of the wavepacket is expected to exhibit 
a relaxation to a steady-state distribution in the limit 
t — > oo. On the basis of numerical calculations it was 
shown in |]l2| that for the time reversal symmetry the 
asymptotic profile f s {n) = \c n (t — * oo)| 2 is given by the 
following expression: 
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where Zoo is the averaged (over energy) localization 
length. Equation (||) was later proved analytically, for 
quasi- ID systems with preserved and broken time rever- 
sal symmetry 0|. It is interesting to note that Eq. (||) 
was derived first for continous ID models with white 
noise potential jl5| with Zoo being the mean free path. 
Remarkably, in spite of the relevant difference between 
the ID and quasi-lD case (the mean free path is of the 
order of b in the latter and thus much smaller than the 
localization length Zoo ~ b 2 ), the asymptotic shape re- 
mains the same in both cases. From (|2j), one finds that 
close to the origin, 

f s (x) ~ exp(-\x\/loo); \x\ < (3) 



while the asymptotic decay is described by 



f s (x) ~ \x\-V 2 expHxl/4/oo); M » 4Zoo, (4) 



revealing that the decay rate s{x) = dlogf s /dx changes 
by a factor 4. 
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FIG. 1. Arithmetic average of the asymptotic wavefunction fs(n). 
The bandwidth is b = 15, and the snapshot is taken at t = 10 4 . 
The data correspond to various crossover parameters: (a) a = 0; 
(b) a = 0.02; (c) a = 0.085; and (d) a = 0.7. 

In Fig. 1 we report some asymptotic profiles f s (x,a) 
for b = 15 and various values of the time-reversal symme- 
try breaking parameter a. These data strongly suggest 
that the decay of f s (n) in the vicinity of the origin is def- 
initely faster than in the tails. As no analytical results 
are available for the crossover regime, it is very tempt- 
ing to compare our numerical results with the theoretical 
dependence (^) derived for quasi-lD systems with pure 
symmetries. The profiles reported in Fig. 2a are the re- 
sult of an arithmetic average over many realizations of 
the asymptotic profile for b — 5, 10, 15 and several a. 
They are plotted with the scaling assumption 



fs(n,t-> oo) = loof s (x); 



n/l c 



(5) 



where Zoo is determined by a fit according to Eq. (g). The 
very good agreement between the numerical results and 
the analytical curve over a broad range of a;- values sug- 
gests that a properly modified theory including the effect 
of an intermediate magnetic field should be able to ac- 
count for the asymptotic profile of wavepackets in quasi- 
lD systems. Furthermore the scaling relation (|J) implies 
that the wavefunction in the crossover regime shows the 
same gross structure (envelope) as in the pure symmetry 
cases, on scales comparable to the localization length. 
This is in contrast to the occurrence of a second local- 
ization scale for the far tails proposed in @J^]. Hence 
we expect that all the scaling laws for the eigenstates 
dominated by the fluctuations of the "envelope" (e.g., 
moments like the inverse participation ratio), also hold 
for the crossover regime. 
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FIG. 2. (a) Asymptotic average profile of the wavepackets for 
b = 5, 10, 15 and various crossover parameters a rescaled according 
to Eq. (^). The light smooth line is the theoretical expression Q). 
(b) The localization lengths lea (open symbols) and (full sym- 
bols) as a function of the time reversal symmetry breaking parame- 
ter a. Various symbols correspond to various values of b; diamonds 
to b = 5, circles to b = 10 and triangles to b = 15. 

In Fig. 2b the asymptotic localization length 1^, as de- 
termined from the best fit with (||) is plotted versus the 
crossover parameter a. One can clearly see that the tran- 
sition from preserved to broken time-reversal symmetry 
is rather smooth and thus we have 



Ma) = (0) 



(0) 



where /3(a) is a smooth function that interpolates be- 
tween the values 1 and 2 for preserved (a — 0) and totally 
broken time reversal symmetry (a ~ 1) respectively. In 
Fig. 2b we also present the localization length derived 
by a direct fit of the average of the logarithm of the 
asymptotic wavefunction — limAr^ 00 (ln(/ s ))/(2iV). 
It clearly shows the same smooth behaviour as Zoo- 

Our numerical results are in contrast to the supersym- 
metry results [[7||| and agree nicely with the transfer ma- 
trix calculations ||. Moreover, they give a definite an- 
swer to the question, whether the assumed discrepancy 
is due to the fact that in the logarithm of the av- 
eraged wavefunction was investigated, while in || it was 
the logarithmic average. Based on our calculations we are 
now able to exclude both possibilities suggested in ||[| ; 
namely that one could attribute the two-scale localization 
phenomenon to anomalously localized states, that domi- 
nate the asymptotic average profile f s (n) or to the non- 
applicability of the Borland conjecture in the crossover 
regime. 
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FIG. 3. Distribution of y = ln/ s (n = 6b 2 ± 2) for b = 5 and 
various values for the time-reversal symmetry breaking parameter: 
(a) a = 0.003; (b) a = 0.03; (c) a = 0.3; (d) Variance of y for 
6 = 5, 10, 15 and various a. The data are rescaled according to 
Eq. (§). 

As mentioned above, the localization of all eigenstates 
implies that for t — > oo the quantum steady state f s (n) 
is localized and fluctuates around the average profile (|^). 
Thus, one can ask about the distribution of f s (n). Hav- 
ing in mind that the asymptotic profile is exponentially 
localized (see Eq. (§)), we computed the distribution of 
the logarithm of the wavefunction y = ln(/ s ) for various 
values of the crossover parameter a. In Figs. 3a-c we re- 
port some representative distributions V(y), which refer 
to the case 6 = 5 with a = 0.003, 0.03, and a = 0.3. In 
all cases, we have found that the distribution of the loga- 
rithm of the wavefunction V(y), is a Gaussian with good 
accuracy. This is in perfect agreement with the analytical 
calculation of Q for V(y). In particular a useful indicator 
is, the spatial growth of the asymptotic variance 



a 2 (n) = ((ln/ s (n)) 2 )-(ln/ s (n)) 5 



(7) 



The results for the cases b — 5, 10 and various val- 
ues < a < 1 are reported in Fig. 3d under the scaling 
assumption 



(8) 



a (n) = a (a;)/v ^oo! n/l 



This further confirms that not only the mean asymptotic 
profile f s (n) but also higher moments change smoothly 
as a function of the crossover parameter a. 

A global characterization of the diffusion and localiza- 
tion properties of a wavepacket is provided by the evolu- 
tion of the mean-square displacement 



m 2 \c rt 



(9) 
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where (•) denotes the average over different realizations of 
the disorder. On the basis of the localization properties of 
the asymptotic wavefunction (see Eqs. (||,||)) one expects 
that, for b ^> 1, M x grows as M x ~ /00(a) 2 . The quan- 
tity M reaches its maximum value at t ~ to- Up 
to that time the evolution of the packet is diffusive, i.e. 
Moo — M(tu) ~ Dtjj, where D is the diffusion constant. 
From supersymmctry |l(| it is known that Zoo A = D 
with A being the mean level spacing. The latter de- 
pends on a as A(a) = VT+ 0.5a 2 A(0). Accordingly the 
following relations hold: 



£^ =/3(a)Vl + 0.5a» (10) 



fop*) _ 0(a) 
t D (0) ~ Vl + 0.5a 2 

where t D (0) ~ 6 3 / 2 and D(0) ~ 6 5/2 @. In the inset of 
Fig. 4 we show the scaling of the diffusion constant D(a) 
for b = 30,40 and various a. Our numerical results are 
in perfect agreement with (|l0|). Equation ( fTo|) suggests 
that the mean-square displacement M(t) follows the scal- 
ing relation 



M(t,a) = M 00 {a)M{t/t D (a)). 



(11) 



The numerical results obtained for the cases 6 = 5, 10, 15 
and various values of the crossover parameter a are re- 
ported in Fig. 4 according to the above ansatz. The 
close coincidence of data reveals the presence of a scaling 
regime already at moderately large 6-values. 
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FIG. 4. Mean-square displacement M(t) for b = 5,10,15 and 
various values of a, confirming the scaling law ([n]). In the inset 
we present the diffusion constant D(pt) as a function of the time- 
reversal symmetry breaking parameter a for b = 30 (squares), 40 
(filled circles). The solid line is the theoretical expectation (|uj)- 

In conclusion, we have studied scaling properties of 
both diffusion and strong localization of wavepackets in 
quasi-lD (i.e., in ID random media with long-range inter- 
actions) . Our numerical data show that the decay rate of 



the asymptotic wavefunction profile changes smoothly, in 
contradiction to the supersymetry results . We have 
also found that the asymptotic shape of the wavepackct 
is well reproduced by the analytical expression (|^) de- 
rived for the pure symmetry cases (i.e. a = 0, 1). More- 
over, in agreement with we find, that the logarithm 
of the asymptotic wavefunction is normally distributed 
for every a, and the variance is scaled according to 
Eq. (||). Another issue addressed in this Letter concerns 
the wavepacket evolution. In particular, we found how 
the diffusion constant D(a) and the diffusion time tr> (a) 
change as a function of the time reversal symmetry break- 
ing parameter a. As a result we were able to establish a 
scaling law ([ll]) for the variance M(t) of the packet. 

We are grateful to Y. Fyodorov and A. Kolesnikov for 
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